The effects of optically induced non-Abelian gauge field in cold atoms 
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We show that N — 1 degenerate dark states can be generated by coupling iV-fold degenerate 
ground states and a common excited state with N laser fields. Interferences between light waves 
with different frequencies can produce laser fields with time-dependent amplitudes, which can induce 
not only U(N) non-Abelian vector fields but also the scalar ones for the adiabatic motion of atoms 
in such laser fields. As an example, a time-periodic gauge potential is produced by applying specific 
laser fields to a tripod system. Some features of the Landau levels and the ground-state phase 
diagram of a rotating Bose- Einstein condensate for a concrete gauge field are also discussed. 



PACS numbers: 32.80.Pj, 32.80.Lg 



5h 



T3 

£ 

O 

> 
00 

o 

o 

-i— > 



i 

-a 
c 

o 
o 



X 



I. INTRODUCTION 

In recent years, cooling atoms in laboratory is open- 
ing up a new playground gathering various ingredients 
and distinct features used to occur in condensed mat- 
ter physics. Most of the prepared systems are described 
by Hamiltonians formally identical to those for electrons, 
which allows one to study the equivalent situations under 
well defined controllable conditions. The electromagnetic 
field, as a kind of gauge field, is known to be playing a 
versatile role there. Like the emergence of the Coriolis 
force in a rotating frame, rotating cold atoms trapped 
in a potential is a useful way to provide an "artificial" 
magnetic field Although it had been noticed two 
decades ago that the Berry phase 0, Q arising in the 
adiabatic dynamics of quantum mechanical systems can 
be regarded as a gauge field, the propagating of slow light 
in a degenerate atomic gas was recently proposed to in- 
duce an effective magnetic field for an electrically neutral 
system |4]. Laser beams with orbital angular momenta 
have been suggested to induce either uniform or radially 
dependent fields [§], which makes it possible to discuss 
the quantum Hall effect and its incompressible nature in 
a degenerate gas of atomic fermions trapped in a two di- 
mensional (2D) confinement [6j. Alternatively, a method 
to realize an artificial magnetic field for neutral atoms in 
a 2D optical lattice was proposed [7[. 

The realization of non-Abelian gauge fields in cold 
atomic systems was theoretically suggested very re- 
cently [TH EH ■ One scheme |T(| employs the laser- 
assisted-state-sensitive tunnelling for atoms with multi- 
internal states, while another one ll| is based on doubly 
degenerate dark states formed via a tripod-coupling sys- 
tem which is a developed scheme of the dark state po- 
lariton As we known, the spatial components of 

gauge potentials were obtained directly whereas the time 
component was introduced through a projection proce- 
dure [lf[ which is normally of higher order. 

The possible emergence of non-Abelian gauge fields in 
various branches of condensed matter physics has ab- 
sorbed much attention recently. For example, Rashba 
and Dresselhaus spin-orbit couplings in certain semicon- 



ductors were perceived to be non-Abelian gauge fields 
[L3| . It is therefore worthwhile to investigate the emer- 
gence of non-Abelian gauge structures in neutral atoms 
systematically. 

In this paper, We show that both the non-Abelian vec- 
tor potentials and scalar ones can be induced in neutral 
atomic systems as long as the amplitudes of the applied 
laser fields are time-dependent. Taking a tripod system 
as an example, we give time-periodic vector and scalar 
gauge potentials using laser beams with particular time- 
dependent amplitudes. We also discuss some concrete ex- 
amples, such as the features of the Landau levels and the 
ground-state phase diagram of a rotating Bosc-Einstein 
condensate (BEC) system in a definite gauge potential. 
We find the gauge potential can make the phase diagram 
of a rotating BEC much richer even if the expression of 
the gauge potential is simple. In next section, we consider 
a tripod system in the laser fields with time-dependent 
amplitudes. In section III, we consider Landau levels and 
the ground-state phase diagram of a rotating system in a 
definite gauge potential. The main conclusions are sum- 
marized in the last section. 



II. MODEL AND GENERAL FORMULATION 

We consider cold atomic ensembles moving in a radi- 
ation field which couples resonantly the atomic ground 
states of iV-fold degeneracy | a) (a = 1, 2, • • • , N) and a 
specific excited state | e) . In terms of the Rabi frequen- 
cies, we have the following Hamiltonian in the represen- 
tation of Hilbert subspace of internal atomic levels 



N 



\a)(e\ 



H.c. 



(1) 



The Rabi frequency depends on the amplitude of the laser 
field where the atom locates. The dressed state reads for- 
mally \ip) = (£ — Hq) 1 !}' \ ip) with Ho the bare atomic 
Hamiltonian, which gives rise to \ip) = J2 a aa I fa) an< ^ 
Wa)=\ a ) + { £ - H y 1 PH' | (p a ) if the operator P pro- 
jecting states outside of the subspace for the degenerate 
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ground states {| a)} is introduced. By iteration proce- 
dure one reaches a formal expression up to any expected 
order: 

[{a I H' | a') + {a\H'(£- H^PH' \ a') + ■ ■ ■] a a , 
= (£ — £o)a a , 

where £o is the unperturbated ground-state energy, i.e., 
Ho \a) — So \a). Since the lowest nonvanishing pertur- 
bation for the Hamiltonian fl} is of the second order, we 
have 



^£ ^2 w* a w a >a a , = £ (1) a a 



(2) 



where A£ is the energy difference between the excited 
state | e) and the degenerate ground states {| a)}, and 
£ = 8q + E' 1 '. By solving this secular equation, we 
get N — 1 degenerate eigenstates | di) with £W = 
and a single non-degenerate eigenstate | s) with £W = 
Y, a h 2 \w a \ 2 /A£ , namely, 



N 

E- 



ii'-' 



•'TV I 



N 



\di) = J2< I")' (* = 1,2,--- ,JV-1). (3) 



The expanding coefficients in the state | di) yield 



w a a n 



0. 



(4) 



they depend on the amplitude of the radiation field. This 
brings about (N — l)-fold degenerate dark states which 
depend on both time and spatial coordinates. 

We apply such laser fields that the Born-Oppcnhcimer 
expansion is applicable. Substituting the state of the 
atom expanded in terms of the dark states (as adiabatic 
basis) 

JV-l 

|¥(r,t)) = \di(r,t)) 

i=l 

into the Schrodinger equation for its translational mo- 
tion, we obtain that 



(p - A) 2 



2M 



^o + V 



(5) 



where ip — (ipi, tp2> " • • > iPn-i) T anc ^ ^ e § au S e potential 
(Ao, A) and V are matrix- valued. Their matrix elements 
are given by 



A,j = ih(di(r,t)\ V \dj(r,t)), 
^o)ij = ih(di(r,t)\ d t \dj(r,t)), 
Va = (di(r,t)\ V(r) \dj(r,t)), 



(6) 



clearly, At = 



A M , (/x = 0,1,2,3). Note that the scalar 
potential in Eq. |G]) differs from the one introduced in 
Ref. [ll| where the scalar potential was obtained by 
means of a projection approach. The Hamiltonian ([5|) 
is invariant under a local U(N — 1) gauge transforma- 
tion. The simplest non-Abelian case is a tripod system 
which provides a U(2) gauge potential. 



which implies that its solutions refer to N — 1 orthogonal 
vectors a = (ai, 02, • • • , a N ) perpendicular to a common 
vector w = [w\, u>2, • • • , w N ). Since (e| H' \ di) = 0, as a 
consequence of (01, those degenerate states | di) are actu- 
ally dark states, which means that the optically induced 
absorbtion and emission occur simultaneously. 

The translational motion of the atom is described 
by the time-dependent Schrodinger equation ihdt 
#(r,t)) = [p 2 /2M + V(r)] \ *(r, t)). Here M denotes the 
mass of the atoms and V(r) an external potential. The 
electric dipole interaction — D-E(r, t) between atoms and 
the laser will not be directly involved in the Schrodinger 
equation although the frequency uj of the laser matches 
the energy difference between the aforementioned ground 
and exited states, i.e., huj = A£. Whereas, we will see 
that radiation fields with time-dependent amplitudes will 
induce non-Abelian gauge fields in the Schrodinger equa- 
tion for translational motions. 

We apply a laser field with bi- frequency oj + and cj_, 
in stead of a single frequency. Here uj = (lu+ + cj_)/2 
matches for the occurrence of transition between the 
atomic ground states and the excited state. The coherent 
superposition of this two waves gives rise to a radiation 
field with a time-dependent amplitude. One can find 
that the Rabi frequencies w a are time-dependent since 



III. EXAMPLE OF A TRIPOD SYSTEM 

As an example, we consider a tripod system which 
provides two degenerate dark states. For convenience, 
we parameterize the Rabi frequencies w a with angle 
and phase variables according to Wi = w sin 9 cos <fie lSl , 
u>2 = w sin 8 sin cf>e lS2 , W3 = w cos 9e 3 , where w = 
| 2 + \w 2 \ 2 + \w 3 \ 2 , tan (f> = \w 2 \/\wi\, and tan 6 = 
\/(\wi\ 2 + |w2| 2 )/|w3| 2 - According to Eq.®, the two 
dark states are given by 

|di) = e~ 4Xl (sin0e~ lSl 1 1) -cos^e"^ 2 |2)), 
\d 2 ) = e^ X2 (cos (9 cos 0e~ 4Sl 1 1) + cos6»sin0e~ iS2 | 2) 
1 3)), (7) 



sin f e 



-iS 3 



where xi an d X2 are arbitrary phase variables. The ma- 
trix elements of the gauge potential (A , A) are calculated 
easily, which yield 

An = h(y(Si + xi) + cos 2 0vs , 2i), 

A12 = /icos6i(isin(20)VS'i2 -iV0)e i(xi " X2) , 

A 22 = ft(cOS 2 9 COS 2 (j)VSi2 + COS 2 0X7523 

+ V(5 3 + X2)), (8) 
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and 



k))n 
k>)i2 
^0)22 



Xi ) 



cos 



hcos6(—sm(2<j))dtSi2 — id t (, 



h(cos 2 9 cos 2 (f>d t S 12 + cos 2 Qd t S 2 3 

+dt(S 3 + X 2)), (9) 



where Sy = Si — Sj. The gauge potentials (Ao,A) cor- 
responding to different \i and X2 are related by a local 
C/(2) gauge transformation. 

We now construct specific radiation fields which lead to 
a time-periodic gauge field. For this purpose, we consider 
two copropagating laser beams in the z-direction with 
slightly different frequencies, whose electric fields read 



E_ 



2 



where k± are the wave vectors of the two fields, u>± the 
frequencies and $± the initial phases. The coherent su- 
perposition of the two fields gives rise to 



E, 



E cos(Afcz - Auit + Ai?)e 



i(k'z+d) -iut- 



where 



2Afc, k+ + fc_ = 2k', u. 



2Au 



lu + +uj_ = 2uj, d + - = 2Ai? and z?+ + = 2d. 
Fixing the initial phase-difference 2Ad on zero, we can 
get a field E stiper with the complex amplitude being 

Ui = E cos( Akz - Aujty {k ' z+§) . (10) 

The Rabi frequencies corresponding to the above field is 

wi = wo cos(Afcz - Auit)e^ k ' z+s) , (11) 

where Hwq = —Eo(D z )i e and (D z )i e = (e \ D z | 1) with 
D z being the z-component of the electric-dipole moment 
of the atom. Additionally, two copropagating and circu- 
larly polarized fields with opposite orbital angular mo- 
menta are applied. The Rabi frequencies of those fields 
are 



w 2 ,3 = w pe 



(12) 



in cylindrical coordinates. 

The gauge potential corresponding to the above fields 
Eq. (fTT|) and Eq.flTJJ) can be calculated from Eq.© and 
Eq.©, 



e ftcos 9 



cos(Afcz — Auit) 
i J cos 2 (Akz - Auit) + p 2 



COS 

cos 8 sin(20) 



cos 9 sin(2</>) 



cos v sin 



sin 



+e z h 



= Hcost 



cosf 



k' sin 2 (f> + k cos 2 <f> 

sin(20)(fc' - k) . Akpsin( Akz - Auit) 
t" i 



i 



2 ' "cos 2 (Akz - Aujt) + p 2 

Auipsm(Akz — Aujt) 



i J cos 2 (Akz - Auit) + p 2 ' 



r 



cosf 



2 "cos 2 (Afcz- Auit) + p 2 

k + cos 2 9 cos 2 4>(k' — k) 



sin(20)(fc' - k) . Akpsin(Akz - Auit) 1 \ 

— 1 

/ 



(13) 



as long as xi an( i X2 are chosen as zero. Here cos <fi = 
cos(Akz — Auit) I \/ p 2 + cos 2 ( Akz — Auit) and cos 9 = 

pj \/2p 2 + cos 2 (Afcz — Auit). We can easily find that the 
gauge potential (A ,A ) is periodic in time with a peri- 
odicity of 2tt I 'Aw. 



tern, the Hamiltonian is written as 



H = + \ Muj2 - x2 + \ Muj ly 2 - nL ( z ) + H ^ ( 14 ) 



IV. EXAMPLE OF A ROTATING SYSTEM 

We now take a two-dimensional rotating cold atomic 
system as an example to investigate the effect of the non- 
Abelian gauge fields. For the rotating cold atomic sys- 



in a rotating frame. Here ui±_ denotes the frequency of the 
trapping potential, the rotating frequency of the sys- 
tem, L(z) = (rxp) z the conventional angular momentum 
operator and Hi nt the interaction between atoms. The 



4 



application of laser fields 



W'l 



w 2 



f hk-2-fx l( _ ky+6i) 

V 2hk 

A(-ky+S 2 ) 

V 2hk 

2 1 X r i(-kz+S 3 ) 

hk 



w 3 = 



with Si, 62 and ^3 being the initial phases and k the 
wave vector of the laser fields, introduces the non-Abelian 
gauge potential A into the above Hamiltonian for rotat- 
ing Bose-Einstein condensates, accordingly 



H = 



(P A) 2 l M 2 2 + i M a 2 

2 - 1 2 A - 



2M 
-O ■ (r x (p 



A)) z + H int . (15) 

Here p = {p x ,p y ), A = (A x .A y ), A x = -\^])t z + \^yl 
and A y = ^■jxt z — ^jxI with f z being the z-component of 
Pauli matrix. Collecting the harmonic trapping potential 
terms into the quadratic term (p — A) 2 /(2M), we can 
write Eq. (fT3j) as 

// JL(p A A)" I U'. ft)L(z) + H mt , (16) 

where A = (A x ,A y ), A x = —Mu^y, A y — Mui ± x and 
L(z) = L(z) - (r x A), = (r x (p - A)),. 

In the limit of O — > ojj_, the Hamiltonian is formally 
the same as that for a charged particle in the background 
of conventional and U (2) magnetic fields. In this case, 
(A x , A y ) becomes (— MQy, MQx), and the correspond- 
ing magnetic field reads B z = 2MQ. The i-component 
of the dynamical momentum is given by 



M 
ih 



H 



= Pi 



A,, - Ai 



(17) 



The commutator between different components of the dy- 
namical momentum is shown to obey 



- X , Ji y\ 



ih(B z +M z ) = ih 



£1 
B 2 



(18) 



where B z = ^t z — 7I refers to the U(2) magnetic field 
which is defined by the corresponding U(2) gauge poten- 
tial. In general, 1* = \eij k {djk k - d k Aj - ^[Aj,A k ]). 
Clearly, different components of the dynamical momen- 
tum are noncommutative as long as either the U(2) mag- 
netic field M z or the conventional magnetic field B z is 
not zero. Following the strategy of algebraic approach to 
Landau levels, we introduce the pseudomomentum whose 
i-component is Kj = f>i + Ai + . The Hamiltonian (in 
the absence of interactions between atoms) and the an- 
gular momentum operator can be written as 



H 



— (P 2 + P 2 ) 
2M [ x+ * h 



1), (19) 



where {B z + M Z )(B Z + B 2 ) _1 = I. We can thus define 
two sets of creation and annihilation operators with the 
above dynamical momentum and the pseudomomentum. 
Each operator has two species 



di 
d 2 



61 
b 2 



= (P x + iPy) 



1 



^/TiBl 




1 



(K X - iKy 







^/2EB 2 


1 



In terms of those creation and annihilation operators, the 
Hamiltonian H = (P 2 + P 2 )/2M is expressed as 



(20) 



whose eigenfunctions and eigenvalues are (| m,mi),\ 
n 2 ,m 2 )) T , and E niU2 = hwi{2ni + 1) + hw 2 (2n 2 + 1), 
respectively, where lj\ — luj_ and lu 2 = uj±_ — 7/M. Here 
the states |ni,mi) and \n2,rri2) are Landau levels with 
n\,m\ and n 2 ,m2 being integers. The integers ri\ and 
TI2 are the Landau level indices, and mi and m 2 label the 
degenerate states within the same Landau level. 

In the following discussion, we assume that the system 
is in the lowest Landau level (LLL), i.e., m, n 2 = 0. The 
assumption is valid in the limit of f2 ~ The wave 
function of the system is written as 



4-1 

* 2 



# = 

*i = E 

mi 

* 2 = E 



mi vmxf. 



m 2 im 2 ip 



(21) 



where a\ = WH/ (Mwi), and a 2 = ^/h/(Mu> 2 ). The 
expectation value of L(z) for the wave function (|2~Tj) is 
given by 



<L(z)) = h /[(r/ar) 2 - l]|*i| 2 dr 
-h I [(r/a 2 ) 2 - 1] |* 2 1 2 dr 
r 2 \^ 2 \ 2 dr. 



(22) 



Then the expectation value of Eq. (fl6|) can be evaluated 
(H) = {lj ± - n)Mcu ± ((r 2 ) 1 + (r 2 > 2 ) + 2MI + H int , 



Hi, 



+Si 2 |*i| 2 |* 2 | 2 



(23) 



5 



where (r 2 ) 1 = /r 2 |\I>i| 2 dr and (r 2 ) 2 = j r 2 \^ 2 \ 2 dr. 
The strengths of the effective two-body interactions are 
.91 = 4nh 2 a s iN/ZM, g 2 = 4nh 2 a s2 N / ZM , and g 12 = 
4irh 2 a s i 2 N j ZM . Here the variation of the trapping po- 
tential along the z-axis is neglected, and the density of 
atoms along z-axis is assumed to be a constant N/Z . 




0.5 



(b) 



FIG. 1: (color online) Ground-state phase diagrams with dif- 
ferent choices of the parameters for a rotating BEC in the 
presence of gauge potentials in the £-Q/uu± plane. The bound- 
aries separating different phases are plotted, cr^, <Ti2, Q/ui± 
and £ are all dimensionless quantities. 

For sufficiently weak interactions, the energy of 
atoms in the rotating frame is minimized when 



the wave functions ^ 2 are composed of 

one component, respectively, i.e., ^ mi ,m 2 



only 



(r 



mi zmiip 



r m 2 im 2 p 



l 2 <) T p 



In this 



case the energy of the system becomes 



(1 



+01 



-) 



mi 



m 2 + 1 



(2m!)! 



+012 



2 2m l(mi !)2 

(mi + m 2 )!(l 



02(1-0 



(2m 2 )! 



2 2m *{m 2 \) 2 

£\(m 2 + l) 



2mi+m2- 



1 m\\m 2 \ 



(24) 



in unit of hujj_. Here £ = j/(Moj±), cr^ = a si N j ' Z and 
(T 12 = a s \ 2 N / Z. For definite values of ft and £, the en- 
ergy of the system depends on the quantum numbers mi 
and m 2 . The optimal values of mi and m 2 corresponding 
to the ground state can be determined by minimizing the 
energy Eq. ([2"3]1. 

From Eq. (|2"4")1 we can see that £i and £o,i become 
lower than £q.o a t critical frequencies of rotation given 



by 



^1,0 
^0,1 



1- y- (I-O012, 

1 _ il-t)*o* _ (25) 



respectively. The vortex begins to be created in one 
component of the system when the frequency of the ro- 
tation reaches the minimal value of {flifl, ^o,i}- For 
some values of £, vortices can be created in both compo- 
nents of the system when Q reaches some critical value 



Thus the state <3> T 



with mi and m 2 being 



both larger than zero becomes more favorable in energy 



for fl > fl r , 



This critical frequency fl r 



calculated from the two conditions. 



^mi,m2 ^mi — 1,7712 



= 0, E n 



C-/7 



can be 



-i=0. 



Thus we can plot out the ground-state phase diagram of 
the rotating BEC in which the boundaries between differ- 
ent phases are obtained by the relation £ mi }Tn2 —£ m ' im ' 2 — 
0. Fig. Q] is plotted for two different choices of the 
parameters where the cr's are small in order to satisfy 
the assumption of weak interactions. One can see that 
the shapes of the phase diagram for the two parameter 
choices change distinctly. This is due to the differences 
of the scattering lengths can break the spatial symmetry 
of the ground state. From Fig.[TJ we can find the angular 
momentum of the system depends not only on the rotat- 
ing frequency ft but also on the strength of the gauge po- 
tential for given interactions. Discontinuous transitions 
may be observed through changing the value of £ even if 
fl is a constant. The presence of gauge potentials makes 
the phase diagram of the system much richer. 

Since we have the expressions of the energy levels near 
the ground state, the phase diagram can be verified with 
the help of the derivative of the ground-state energy. In 
Fig. [21 we plot the derivative of the ground-state energy 
along the (1, —1) direction versus (£, (l/w±), 



d£ 



d£ a 



d£ d(Sl/u 



(26) 



where £ g is the ground-state energy of the system. One 
can see that discontinuous steps emerge in the F'(£, O) 
surface (i.e., singularities for the derivative of F'(^,fl)). 
This exhibits the boundaries separating different phases 
as it is also an effective quantity to characterize quantum 
phase transitions [l5| . These boundaries indeed concise 
with those given in Fig. [TJ As well known, there are var- 
ious quantities for characterizing quantum phase tran- 
sitions, such as the quantum entanglement [16| and the 
fidelity [l7| . The former is convenient if the reduced den- 
sity matrix for paring correlation is computable while the 
latter is convenient if the ground state is expressible. In 
present case, the derivative of the ground-state energy is 
the most convenient quantity for that purpose. 

The above discussion is based on the assumption of 
weak interactions, which requires that the interaction en- 
ergy must be relatively much smaller than kinetic and 
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FIG. 2: (color online) The curve surfaces are the derivatives 
(in the unit of hui±) of the ground-state energy along (1, —1) 
direction versus (£,Q/uj±) for parameters a\ = 0.2, 02 = 
0.4, CT12 = 0.3 (top panel) and o\ = cr 2 = (T12 = 0.4 (bot- 
tom panel) respectively. The discontinuous steps precisely 
occur at boundaries separating different phases illustrated in 

Fig.m 

trapping potential energies for each component. The for- 
mer is of order ~ pigi (i=l,2), and the latter is of order 
~ hu)i, where pi is the particle density of the zth com- 
ponent in the x-y plane. The cross section of the ith 



component is iraf in the weak interaction limit. Thus, 
the density pi is l/naf. So we can obtain the effective 
coupling constant <y% <C 1 from the condition pigt <C huii. 

If the value of Oi violates the condition of weak inter- 
actions, the energy of the system will be minimized when 
many components in the expansion of Eq. (|21|) contribute 
to the wave function ^j. In this case, the equations (j2"3")l 
is still valid. To calculate the expectation value of r? 
and interaction Hamiltonian iJ, nt given in Eq. (f2"5)) . the 
numerical simulation and the averaged vortex approxi- 
mation [18[ are efficient methods. Note that two sets of 
vortex lattices may be formed in this case. The structure 
of vortex lattices depends on the strength of interactions 
as well as that of the gauge potentials. The presence 
of gauge potentials can increase the diversity of vortex- 
lattice structures, so it is worthwhile to further study 
these systems. 



V. CONCLUSION 

In summary, we have given a general description of the 
dark states and found that N — 1 degenerate dark states 
can be generated with the help of N laser fields cou- 
pling ./V-fold degenerate ground states with a common 
exited state. Interferences between two waves with dif- 
ferent frequencies can produce such radiation fields that 
their amplitudes are not only spatially dependent but 
also time-dependent. They can induce vector and scalar 
gauge potentials. As an example, we have considered a 
tripod system for which one can obtain a time-periodic 
gauge potential using a specific laser field. We have given 
a configuration of laser field that leads to a uniform U(2) 
magnetic field, in which we have discussed the features of 
Landau levels and studied the quantum phase transitions 
of a rotating BEC using the derivative of the ground-state 
energy. We have shown that the presence of gauge poten- 
tials will make the picture of rotating BECs to be much 
richer. 

This work is supported by NSFC No. 10225419 and 
No. 10674117. 
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